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ABSTRACT 

We present circular and linear cross-correlation tests and the "friend-of-friend" anal- 
ysis for phases of the Internal Linear Combination Map (ILC) and the WMAP fore- 
grounds for all K-W frequency bands at the range of multipoles I < 100. We compare 
also Tegmark, de Oliveira-Costa and Hamilton (2003) and Naselsky et al. (2003) 
cleaned maps with corresponding foregrounds. We have found significant deviations 
from the expected Poissonian statistics for all the cleaned maps and foregrounds. Our 
analysis shows that, for a low multipole range of the cleaned maps, power spectra 
contains some of the foregrounds residuals mainly from the W band. 

Key words: cosmology: cosmic microwave background — cosmology: observations — 
methods: data analysis 



1 INTRODUCTION 

The recently-published Wilkinson Microwave Anisotropy 
Probe ( WMAP) data sets (see Bennett et al. 2003 a-c, Hin- 
shaw et al. 2003 a-b) open a new epoch for the CMB in- 
vestigation. These data serve for the development of more 
refined technique for future high resolution measurements 
and for a choice of the realistic cosmological model. 

The WMAP has observed the full sky in five fre- 
quency bands: K (centered frequency 22.8 GHz), Ka (33.0 
GHz),Q (40.7 GHz), V (60.8 GHz), and W (93.5 GHz) and 
produced five maps represented in n = 12 x 512 2 HEALPix 
(Gorski et al. 1999) pixels. Then, using smoothing of the 
maps to 1° and performing cleaning of the combined mul- 
tifrequency map by minimization of the rms variance for 
each pixel, the WMAP team produced the Internal Linear 
Combination (ILC) map taking 12 optimization coefficients 
Wi,i — 1,2,..., 5 for 12 disjoint sky regions. Thus, the ILC 
map has a minimal variance corresponding to the minimiza- 
tion of the Galaxy and foreground contamination at the 
range of multipoles £ < 100. 

Tegmark, de Oliveira-Costa and Hamilton (2003) (here- 
after TOH) have suggested to use ideologically similar but 
technically different optimization scheme based on Tegmark 
and Efstathiou (1996) (hereafter TE) method. The main 
goal of the TOH method is to combine the K-W band maps 
into the Foreground Cleaned Map (FCM) using variable 
weighting coefficients u)i(l) for each frequency band i=K- 
W. In order to minimize the total unweighted power from 
the Galaxy, foregrounds and noise separately for each har- 



monic £, TOH subdivided the sky into 9 disjoint regions 
and performed cleaning of the 5 frequency maps in each 
region. Unlike the ILC map, the TOH FCM has smaller 
power at the multipole range £ < 12. Both the ILC and 
TOH FCM are available from the web and can be used 
for analysis of the foregrounds and the statistical properties 
of the cleaned maps. For example, Komatsu et al. (2003) 
tested the non-Gaussianity of the WMAP CMB signal us- 
ing Minkowski functionals as a statistic. Chiang et al. (2003) 
have also tested the statistics of the phases for the ILC and 
the TOH FCM maps and discussed non-Gaussianity of the 
TOH FCM map. Dineen and Coles (2003) have presented 
a diagnostic of the Galactic synchrotron contamination and 
have discovered the cross-correlation between Faraday rota- 
tion measure and the TOH FCM map (see also Coles et al. 
2003). 

Recently Naselsky et al. (2003) proposed the Phase 
Cleaning Method (PCM) for decomposition of the CMB 
signal and foregrounds using the WMAP data. The main 
idea of the PCM is to minimize the variance of the com- 
bined K-W maps and cross-correlations of the CMB and 
foreground phases simultaneously. The result of PCM in its 
application to the WMAP data reproduces well the power 
spectrum of the best-fit WMAP ACDM model at the range 
£ < 50. The main target of our paper is to extend the res- 
olution of the derived PCM CMB map up to £ ~ 100, to 
perform the the cross-correlation analysis (circular and lin- 
ear statistics and "friend-of-friend" analysis) of the CMB 
and foreground phases and to estimate their interconnection. 



2 Naselsky, Doroshkevich & Verkhodanov 




Figure 1. The map for reconstructed CMB at the range of multipoles i < 100 (top left) and differences between ILC and 
PCM maps (top right), TOH FCM and PCM (bottom left) and TOH Wiener and PCM maps (bottom right). 



We have included the ILC , TOH FCM and the PCM maps 
in our analysis and have discovered significant deviation of 
the phase statistics from the expected for uniformly random 
phases. We have shown that such a kind of non-Gaussianity 
arises owing to the contamination of the foregrounds signal 
in the ILC, TOH FCM and the PCM maps. 



2 PHASE CLEANING METHOD AND THE 
CMB EXTRACTION 

The fluctuations of measured CMB plus foregrounds radia- 
tion on a sky sphere can be expressed as a sum over spherical 
harmonics: 

oo t 

AT(6, <p) = T(6, <p)-(X)=^2 a e™Y lm {6, <p), (1) 

£=0 m=-l 

where ae m are the coefficients of expansion, (T)=2.73K, 
and (AT)=0. Homogeneous and isotropic CMB Gaussian 
random fields (GRFs), as a result of the simplest infla- 
tion paradigm, possess ai m modes whose real and imaginary 
parts are independently distributed. The statistical proper- 
ties of this field are completely specified by its angular power 
spectrum Ct, 
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which are uniformly distributed at the range 0, 27r. For the 
foregrounds, the signal is obviously non- Gaussian. So, for 
the combined CMB + foregrounds signal we define 



where \ae m \ ls the modulus and \&£ m is the phase of each t, m 
harmonic. In practice, each of the WMAP K-W maps is de- 
composed into set of a-* where a = ae m an d index j — 1, 2. .5 
corresponds to K, Ka, W bands, using the HEALPix code 
(Gorski et al. 1999). 1 

The basic idea of the PCM is to generalize the TOH and 
TE minimization scheme including also the minimization 
of the cross-correlations between derived CMB signal and 
foregrounds (Naselsky et al. 2003). The method does not 
require any galactic cut-offs or disjoint regions. Following 
Naselsky et al. (2003), we consider the combinations of the 
WMAP maps: Ka-Q, Ka-V and Q-V, for which the higher 
correlations between foreground phases takes place. 2 . For 
each pairs of the maps, we introduce weighting coefficients 
w^(^) similar to TOH and TE methods and minimize the 
variance of the derived map per each mode t. Neglecting 
the instrumental for I < 100 noise and taking into account 
the beam shape deconvolution, we consider the signal as a 
superposition of the CMB and the foregrounds. 

From all a^\ we have found the phases = 
(sec Eq ©)• Each set of the a^' coefficients is defined 
by a combination of the different foreground coefficients 

i2 k a k ) = G<J) and the CMB si s nal &cmb ' where index k 

marks the synchrotron, free-free and dust emission. For each 
combination of the maps, the derived CMB map is defined 
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In addition we used the GLESP pixelization scheme (Doroshke- 
vich et al. 2003a) in order to control the accuracy of the a.U) 
harmonics estimation. 

2 We did not include W and K bands to the separation procedure 
because of the peculiar phases (Naselsky et al. 2003). 
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Figure 2. The power spectrum for the PCM reconstructed CMB 
(thick solid line), ILC (dash line), TOH FCM (dash dot line) and 
TOH Wiener filtered map (thin solid line). 
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and 3"' is the phase of j'-channel. 

The phases of the reconstructed CMB are related to the 
foregrounds amplitudes and phases as follows (Naselsky et 
al. 2003) 

V .w (3, |G (j) |sin(* (j) -£)cos(* (M) ) 

VJ> M — t | arcGin 3 (7) 

J2j wO')|GW| cos*0) + |a cm6 | cos£ 

where £ is the true CMB phase and are the foreground 
phases. As is mentioned in Naselsky et al. (2003), the phases 
of the reconstructed CMB signal Of would have correla- 
tions with those of foregrounds. The main goal of the PCM 
is to minimize such correlations by minimizing of the weight- 
ing variance 
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using Eq.JTJ and weighting coefficients w 7 in a form 
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Note that minimization of the weighting variance Eq.© is 
equivalent to minimization of the error bars of the CMB 
reconstruction ^ m \af^ — a,i™ b \ 2 — > min, where a^™ 6 is the 
true CMB (Naselsky et al. 2003). 

Practical implementation of the PCM includes the it- 
eration scheme in which the first step corresponds to min- 
imization of the variance of the derived map in each mode 
£ with the choice of coefficients Eq.l|S|. This step of opti- 
mization reconstructs the pre-CMB and pre-foreground map 



Figure 3. The maps for differences S — ILC — Foregrounds for 
V (top) and W (bottom) WMAP bands. 



for each pair of channels where the foreground maps are a 
simple subtraction ag m coefficients of the K-W signals and 
derived pre-CMB ai m . Then, using the moduli and phases 
of the pre-foregrounds we perform the next iteration, using 
Eqs.©, JSJ and so on. This iteration scheme is stable and 
reproduce quite well the CMB signal from each pair of K-W 
maps after two steps of iteration. 

However, the CMB maps reconstructed from each pair 
of K-W maps are slightly different because of residues of 
the foregrounds. To minimize such residues, as the last step 
of the PCM we use the so called MIN-MAX filter for each 
pixels of the map. The MIN-MAX filter compares signals, 
AT P , in each pixel p in a set of maps. In this method, we 
chose the minimal amplitude of AT V in the pixel as the 
CMB signal, i.e. L min ({AT p W }) -» AT p min , where |AT™ in | = 
min{| ATp^ |}. To estimate residues from point sources, we 
find the maximal amplitude of the signal in each pixel 



L max ({AT p W }) -» AT™ ax , where \AT£ 



maxflAT^I} 



for all the pre-CMB maps produced by the method Eq.@. 
For true CMB, the difference AT p max -AT™ in is equal to zero 
and we consider AT™ ax — AT™ ln as a measure of deviation 
from the CMB map. 

The reason for such a filter is quite obvious. The sig- 
nal in each pixel is a superposition of the CMB signal and 
(small) residues of the foregrounds. If correlations between 
the CMB signal and foregrounds are minimal then we ex- 
pect that all the deviations of the ATp in the pixel are 
caused by the residues. In our case, two approximations of 
the CMB map restored by Ka-V and Q-V channels (Nasel- 
sky et al. 2003) have been used and AT mm has been taken 
as the CMB map. 

The result of the PCM application to the WMAP data 
at the range £ < 100 is shown in Fig. 1 for the CMB signal. 
In this Figure, the differences between the ILC and PCM 
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Figure 4. The circular correlation between the ILC (thick solid 
lines), FCM (thin solid lines) and PCM (dashed lines) cleaned 
signals and derived foregrounds for K— W channels. 



maps (top right), TOH FCM and PCM (bottom left) and 
between TOH Wiener and PCM maps (bottom right) are 
also shown. As is shown in Fig. 1, the main difference is con- 
nected with sources in the Galactic plane and the harmonic 
which is probably related to the galactic emission. The 
difference shows foreground components being absent in the 
PCM map. 

For the simplest estimator of the power spectrum, 
C{£) = (2£ + l)' 1 J2 m \aim\ 2 , we show in Fig. 2 the power 
spectrum for the ILC,"fCM, PCM and TOH Wiener filtered 
maps. As one can see from Fig. 2 the PCM map reproduces 
well the spectrum for the TOH Wiener filtered map. 



3 CORRELATIONS OF THE CLEANED 

SIGNALS AND FOREGROUNDS FOR THE 
ILC, FCM AND PCM MAPS 

Three cleaned maps (ILC, FCM and PCM) are produced 
from different methods of the foreground component separa- 
tion. For the perfect separation of the CMB signal and fore- 
grounds we can expect small random correlations between 
their phases. However, these methods provide the approx- 
imate separation only. Therefore, these remaining cross- 
correlations between phases of the cleaned signal, <j} 3 , and 
the foregrounds, tpf, can be used to characterize the degree 
of separation achieved. 

To do this, here we will use the simplest phase- 
phase circular (Fisher, 1993) correlation coefficients. We 
will use also the linear correlation coefficients and more re- 
fined "friend-of-friend" statistics discussed, for example, by 
Roeder (1992) neglecting circularity of the phases. Indeed, 
equivalence of the phases (f> s &z<f> s ± 2tt and ipf kipf ± 2-7T, 
changes their correlation functions and cumulants and the 
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Figure 5. The circular correlation between the ILC (thick solid 
lines), FCM (thin solid lines) and PCM (dashed lines) cleaned 
signals and free-free emission for K-W channels. 



results obtained with the "friend-of-friend" statistics. To 
suppress this effect we compared our result with 200 ran- 
dom realizations of Poissonian process prepared in the same 
manner (see Sec. 3.3). 

For five frequency channels K-W, the maps are taken 
from the WMAP web site 3 , and all the phases are ob- 
tained by the spherical harmonics decomposition using the 
HEALPix (Gorski et al. 1999) and GLESP (Doroshkevich 
et al. 2003a) codes. We consider separately two ranges of 
multipoles, 2 < £ < 50 and 51 < t < 100. For the ILC, 
FCM and PCM CMB maps, the derived foregrounds are de- 
termined as differences between the signal (S) and cleaned 
map (C) F = S—C for each (£, m) . In addition, for the ILC 
cleaned map, WMAP own foregrounds given in the same 
web-site were also used. Examples of difference of S, ILC 
and WMAP own foregrounds for V and W WMAP bands 
are shown on Fig. 3. 

3.1 Circular cross-correlation of the phases. 

Following Fisher(1993) we define the statistics 



M sp = — !— V] exp[ip((j> m - {4>m))\ , 



M fp = 



t V] exp[ip(ip m - {ip m })] , 

771 

/1-mo.x \ 
{4>m} =tan _1 ^2 sin 0m/ ^2 COS ^m ! 

\ TTl 777, / 

3 http: //lambda. gsfc .nasa.gov/product/map/m_products . cfm 
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Figure 6. The circular correlation between the ILC (thick solid 
lines), FCM (thin solid lines) and PCM (dashed lines) cleaned 
signals and synchrotron emission for K-W channels. 



Figure 7. The circular correlation between the ILC (thick solid 
lines), FCM (thin solid lines) and PCM (dashed lines) cleaned 
signals and dust emission for K-W channels. 
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Table 1. Circular cross-correlation coefficients, (R 3 f), between 
phases of the ILC signals and own and derived foregrounds, and 
PCM and FCM signals with their derived foregrounds. 
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where M sp and Mf p are the p-th trigonometric moments 
of the samples, (<f>) and (tp) are corresponding mean direc- 
tions, R a f{£) is the circular cross-correlation coefficient in 
each mode £ and r s f is the mean circular cross-correlation 
coefficient for all phases. For m = and for all I phases 
4>(£,0) = ip(£,0) — and here we neglect them. 

For the K-W bands the circular coefficients, R s f(£), are 
plotted in Fig. 4 for the ILC cleaned signal and its own and 
derived foregrounds and for PCM signal and derived fore- 
grounds. For the first three channels, these coefficients are 
quite moderate and do not exceed the random scatter (la) 
obtained from 200 random realizations. Note that, for all 



the bands, the shape of the functions R a f(£) are quite sim- 
ilar to each other what reflects strong correlation of phases 
in all the foregrounds (Naselsky et al. 2003) . As one can see 
from Fig. 4, for the V & W channels the cross-correlations of 
both the PCM and ILC phases with the derived foregrounds 
seems to be quite significant. 

The same tendency is seen from the estimations of the 
mean coefficients, {R s f}, for 2 < £ < 50 listed in Table 1 for 
all three cleaned signal and the own WMAP and three de- 
rived foregrounds. For first three channels these coefficients 
are small but they become significant for channels V and W 
and derived foregrounds. 

Correlations of phases of different foregrounds with the 
ILC, FCM and PCM cleaned signals are plotted in Figs. 5-7. 
For the free-free foregrounds, correlations are quite moder- 
ate and they exceed the random scatter only for the 35 - 
36 harmonics. However, for the synchrotron and dust fore- 
grounds we see the significant (95%) correlations for the 21 - 
23 harmonics for all cleaned signal maps. In all the cleaned 
maps this range corresponds to minima of the power (see 
Fig.2). 



3.2 Linear cross-correlation of the phases 

For £min < £ < £max, the linear correlation coefficient be- 
tween phases of the cleaned signal, </> B , and foregrounds xpf, 
is defined as follows: 



Tsf (£min , £max ) 



(11) 
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0.088 
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0.033 
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Table 2. Linear correlation coefficients between the phases of the foregrounds and cleaned signal for all the K-W 
bands. The first three columns represent r s f for WMAP own foregrounds and the ILC, FCM and PCM cleaned sig- 
nals. In other columns the cross-correlation coefficient for ILC, FCM and PCM signals with their derived foregrounds 
are listed for two ranges of multipoles, I < 50 and 51 < I < 100. 
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For the foreground phases, tpf, the mean value, (tpf), and 
the variance, <x^, are denned with similar relations. Here Nt 
is the number of phases in the sample under consideration. 
For m = and for all £, phases (j> B (£, 0) = ipf (I, 0) = and 
here we neglect them. 

The linear correlation coefficients for 2 < £ < 50 and 
51 < I < 100 are listed in Table 2 for the ILC, FCM and 
PCM cleaned signals and five frequency channels. First three 
columns characterize the correlation with the ILC own fore- 
grounds. Six other columns represent these coefficients for 
the same signals and derived foregrounds. For six samples 
the correlations between the signal and derived foregrounds, 
r a f(imin,£max), increase with the channel frequency and are 
comparable for a given frequency. The random scatter of 
the coefficient determined by averaging of 200 random re- 
alizations is r s f(£ m in, (-max) ~ 0.06 at 68% CL. As is seen 
from Table 1, for high frequency channels V and W the mea- 
sured r s f(£ m i n , iraax) exceeds this value. This fact indicates 
the noticeable correlation between the cleaned signal and 
derived foregrounds in channels V and W and the limited 
precision of separation methods. 

For all the channels and WMAP own foregrounds, the 
coefficients r s f(£min, (.max) are less than the random value 
that indicates high efficiency of corrections of own fore- 
grounds. However, this result depends upon the foregrounds 
rather than the cleaned signal. For the FCM and PCM 
cleaned signal and the same foregrounds, the coefficients, 
r s f{£min,(max) are also listed in the same Table. For all 
samples these coefficients r s f(£min, (max) do not exceed the 
random value. Therefore, the small correlations between the 
own foregrounds of ILC and the cleaned signals is deter- 
mined by the properties of foregrounds. 

By definition, this coefficient is more sensitive to (max 
because the majority of phases used comes from £ ~ (max- 
However, it decreases when we correlate phases with differ- 
ent I 



3.3 Linear correlation coefficient of the phases 
per each multipole £ 

Linear correlation coefficient of phases of cleaned signal and 
foregrounds, r s f((), can also be found for each £ in the same 
manner as it was been done for the circular cross- correla- 
tion coefficient in Sec. 3.1. In the case, the coefficient is de- 
fined by the same relations 1 11 1 neglecting the summation 
over £ and for Ne = (■ It allows one to characterize prop- 
erties of each mode separately and to determine harmonics 
with maximal r s f(£). However, the small sample statistics 
of phases at moderate £ increases significantly its random 
scatter. Note also that, as is seen from the definition, this 
coefficient characterizes the sample used. Thus, in contrast 
with the circular cross-correlation coefficient, the direct av- 
eraging r 3 f (£) over £ does not reproduce results obtained in 
previous Section. 

The coefficients r s f(t) are plotted in Fig.|H]for channels 
V and W for two foregrounds of the ILC cleaned signal and 
for the PCM cleaned signal. For each £ the error bars are 
found for 200 random realizations and are well consistent 
with the expected behavior ~ ^ _1//2 . 

For channels K-V and for all three samples, the cor- 
relation coefficients do not exceed the random scatter but 
the disposition of their maxima and minima are very simi- 
lar. For the channel W we see significant differences between 
functions r 3 f{l) for different samples. More pronounced cor- 
relations are seen for 4 < I < 10, £ = 13 - 15, £ = 18 - 27, 
( — 32 — 36 and ( = 45 — 50 As is seen from the compari- 
son with Fig.|21 some of these harmonics correspond to local 
minima in Ci. This fact indicates that, perhaps, these pecu- 
liarities are caused by the same sources as the Galactic dust 
emission. 



3.4 "Friend-of-Friend" statistics 

To compare the foregrounds and the foreground-cleaned sig- 
nal with random Poissonian process we use a more refined 
technique, the "friend-of-friend" approach for the 2D phase 
diagramme. In these diagrammes each point has coordinates 
Xi = 4> s ((,m),yi = ipf(l,m) where <p a (£,m) and ipf(£,m) 
are the phases of the cleaned signal and the foregrounds, 
respectively. 

For each diagramme, we find the fraction of clusters, 
fit, with a given richness, M, in wide range of linking 
lengths, ri n k, which determines the maximal point separa- 
tion within a cluster. Here we use the normalized linking 
length, 

link = irr? nk {n), (n) = N phs /4n 2 , (12) 
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Figure 8. Linear correlation coefficient, r a f, between the CMB 
and foregrounds vs. the harmonics number £, in V and W channels 
for ILC (solid line), ILC and its own foregrounds (dash line, which 
overlapped practically the solid line), and PCM filtered map (long 
dashed line). 

where ri n k is the actual distance between two points, and 
Nphs is the number of points used. As was discussed in White 
(1979) and Borgani (1996), the fractions /m depend on the 
correlation functions (or cumulants) of all orders that pro- 
vides high sensitivity to deviations of the sample from the 
Poissonian-like one. This sensitivity is weak for M = 1 & 2 
and fast increases for larger M. But the random scatter of 
/m increases for larger M as well, and the analysis becomes 
ineffective for M > 7. 

For the truly Poissonian sample of points, the corre- 
lations functions of all orders are equal to zero and these 
fractions can be approximately described as follows: 

fM{knk) ~ [1 — exp( — ll„k)] M 6Xp( — ll nk ) , 

where the first term describes the probability to find M 
separation between points < link while the last term gives 
probability to find larger separation for the M + 1 point. 
The scatter of /m is determined by averaging of /m over 
200 random realizations of diagrammes analyzed in the same 
manner. 

To find clusters with a given richness we use the Min- 
imal Spanning Tree technique which, firstly, allows one to 
connect all points within a unique tree. Further on, rejection 
of edges of a tree larger than a chosen linking length trans- 
forms the tree to the system of clusters of all richness for a 
given linking length. The Minimal Spanning Tree technique 
was proposed in Barrow et al. (1985) and van de Weygaert 
(1991) and was applied for analysis of rich galaxy catalogues 
(see, e.g., Doroshkevich et al. 2001, 2003b). 

Firstly, we are testing the hypothesis that the distribu- 
tion of phases are Poissonian-like and correlations are negli- 
gible for both foregrounds and cleaned signal. Main results of 
our analysis are presented in Figs. 9,10 where for all five fre- 
quency channels variations of fractions fM{hnk), 1 < M < 7 



are plotted versus the linking length, k n k- As is expected, the 
functions /m {link) are quite similar to Poissonian ones for 
M = 1 and M = 2, and even for M — 3 the deviations are 
moderate. However, for M > 4 the deviations from the Pois- 
sonian samples become significant, especially for V and W 
channels. This result illustrates the expected non-Gaussian 
character of the cleaned signal and its correlations with the 
foregrounds. 

This result is expectable because of non-Gaussian char- 
acter of the foregrounds signals. However, comparing the 
phase diagramme for the cleaned signal and foregrounds 
with 100 diagrammes prepared for the foreground and ran- 
dom phases we can quantify their divergences using \ 2 
statistics. These results are shown in Table 3. If for M = 1 
and M = 2 we see small \ 2 < N where N is the number 
of linking lengths used then for M > 3, we have \ 2 > N . 
Larger ratio \ 2 /N is found again for higher frequency chan- 
nels V and W. This fact indicates that all the cleaned maps 
contain some residues from the dust emission. 



4 CONCLUSION 

In this paper we have proposed new methods of investiga- 
tion the statistical properties of the signal derived from the 
WMAP K-W bands foreground cleaned maps. With this ap- 
proach we show significant correlations between the phases 
of the cleaned maps and the foregrounds. We have compared 
three cleaned maps, ILC, TOH FCM and PCM, and have 
investigated variations of the linear cross-correlation coeffi- 
cients versus the multipole index, I. We suggest also more 
sophisticated "friend-of-friend" statistics. All the methods 
described in Sec. 3 show significant cross-correlations be- 
tween the phases of the cleaned maps and the foregrounds 
which manifest themselves more clearly for the W band of 
the WMAP. We have pointed out that some of peculiar- 
ities of the ILC and PCM power spectrum, for example, 
local minima at I ~ 6, and t ~ 20 — 25 are accompanied 
by significant cross-correlations between phases of the fore- 
grounds and cleaned signal. Such correlations are important 
indicators for investigation and detection of possible non- 
Gaussianity of the CMB signal which could be a byproduct 
of the component separation methods. 
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Figure 9. Fraction of clusters with richness M = 2 — 7 vs. dimcnsionlcss linking length, for the PCM reconstructed 
CMB (thick solid line), ILC (dash line), TOH FCM (thin solid line). The dots represent a set of random samples. The range 
of multipoles for the phases is 2 < t < 50. 
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Table 3. Comparison of degrees of clusterization between the foreground and random phases (200 realization) and the 
foreground and CMB phases for five channels of the ILC, FCM and PCM samples. The \ 2 an d the number of linking lengths 
used are presented for clusters with richness between 1 and 7. 
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5 APPENDIX 

In this section we would like to introduce an analytical ap- 
proach for investigation of the cross-correlation coefficient 
r s f(£) between foregrounds and derived CMB phases per 
each mode i. We define r s f(l) as 



r k s f (e) 



— Y 



{At) 



where k marks the foreground derived from the WMAP 
frequency band. Here we assume that the distributions of 
both phases are close to homogeneous ones with (^ r )=0, 
(* 2 ) w tt 2 /3 and we average a correlation coefficient over 
m for a given I. For the phase of the cleaned CMB signal we 
will use Eq.lO, which reflects directly the cross-correlation 
between the derived CMB and foregrounds phases. Then, 
using the definition ^ = aj^ k , where is the cross- 
correlation coefficient between the foreground phases in j'-th 
and fc-th channels, from Eq.(Al) we obtain 



'<■/ 



(£) = 



EE- 



(a 



(A2) 



jki 



Let us assume that cross-correlation between different 
foreground phases is very strong and corresponding coeffi- 
cients are aifc ~ l 5 . Then, from (A2) we obtain 



k 

r sf 



{A3) 



For the PCM method, w"' coefficients are the functions of 
£, but not m, while ratio |G W3 !/!a cmb | depends on m for a 
given value £. 

From Eq.JSJ we get 



I Af cmbi 

a - a 



(AA) 



and, because the foregrounds seem to be highly corre- 



lated, 



sign 



X> Cf) |G«|W 

\ j ) 3 



w (i) |G (j) | 



see Naselsky et al. (2003) for details of the phase correlations. 
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q . I Q ( M ) _ -cmil 

*M = -$r t £ ' | a c^ ' ^(|a^>|-|a^|) (A5) 

m— 1 

Thus, the cross-correlation coefficient r^(£) simply de- 
scribes the accuracy of the reconstruction of the CMB signal. 
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